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Abstract — We analyze the practical performance of quantum 
polar codes, by computing rigorous bounds on block error 
probability and by numerically simulating them. We evaluate our 
bounds for quantum erasure channels with coding block lengths 
between 2 10 and 2 20 , and we report the results of simulations 
for quantum erasure channels, quantum depolarizing channels, 
and "BB84" channels with coding block lengths up to N = 1024. 
For quantum erasure channels, we observe that high quantum 
data rates can be achieved for block error rates P e < 10 -4 
and that somewhat lower quantum data rates can be achieved 
for quantum depolarizing and BB84 channels. Our results here 
also serve as bounds for and simulations of private classical data 
transmission over these channels, essentially due to Renes' duality 
bounds for privacy amplification and classical data transmission 
of complementary observables. Future work might be able to 
improve upon our numerical results for quantum depolarizing 
and BB84 channels by employing a polar coding rule other than 
the heuristic used here. 

Arikan's polar codes are an important recent development 
in coding theory (T). These codes provably achieve the Shan- 
non capacity of symmetric binary-input memoryless channels 
and have an O(TVlogTV) complexity for both encoding and 
decoding, where N is the number of channel uses. Polar codes 
exploit the channel polarization effect, in which a particular 
recursive encoding induces a set of virtual channels, such 
that some of the virtual channels are near perfect for data 
transmission and the others are near useless for this task. The 
fraction of the virtual channels that are near perfect is equal 
to the symmetric capacity of the original channel. 

Several authors have now extended the ideas of Arikan 
to the domain of quantum information theory, in order to 
accomplish a variety of information processing tasks, includ- 
ing classical data transmission J2], (5), private classical data 
transmission |4|, (5), and quantum data transmission [6|, |4j, 
0. Among these works, Ref. J6] demonstrated that all of 
the above impressive features of classical polar codes are 
preserved when constructing quantum polar codes for sending 
quantum data over quantum Pauli channels or quantum erasure 
channels. Refs. J4], Q then followed up with similar results 
for private classical data transmission over these channels. The 
importance of Pauli channels stems from the fact that any 
noisy qubit channel can be "twirled" to a Pauli channel with 
a quantum capacity not higher than the original channel. A 
quantum depolarizing channel is a Pauli channel that is often 
used as a "worst-case scenario" noise model, and a quantum 
erasure channel is a simplified model of photon loss. 



In light of the above theoretical works, it seems natural 
now to assess the practical performance of finite-blocklength 
quantum polar codes on any of the aforementioned channels. 
One could get a better idea of this by deriving rigorous 
bounds on performance like Arikan's bounds for classical 
erasure channels 0], and another way to do so is simply by 
simulating their performance numerically. Indeed, since the 
quantum decoder for Pauli and erasure channels is a coherent 
version of Arikan's successive cancellation decoder combined 
with other efficient operations J6], it is possible to have an 
efficient numerical simulation of these codes' performance. 

In this paper, we derive bounds on the performance of 
quantum polar codes on quantum erasure channels, and we 
report the results of numerical simulations determining the 
performance of quantum polar codes on quantum erasure, de- 
polarizing, and BB84 channels. In our numerical simulations, 
we consider the performance of quantum successive cancella- 
tion decoding for block lengths of 64, 256, and 1024. We find 
that with these moderate block lengths, relatively large rates 
can be achieved on the quantum erasure channel with a block 
error rate P e < 10 -4 . For the depolarizing and BB84 channels, 
we compare our results against the well-known hashing bound 
for quantum data transmission. For these channels, we discuss 
how our results serve equally as simulations of polar codes 
for private classical data transmission. 

We structure this paper as follows. To start, Sections [IJ [II] 
and III review some of the recent advances in polar coding 



for quantum channels (2), @, |4j, Q, 0, 0, focusing on 
the approaches given in Refs. 0, 0, 0. Section IV derives 



our bounds on quantum polar code performance for quantum 
erasure channels. We then detail the results of our numerical 
simulations in Section \V\ In Section fVl] we discuss how our 
results apply in the setting of private classical communication, 



and we finally conclude in Section VII 



I. Polar codes for classical communication 

Consider a channel W with a classical input x € {0, 1} and 
a quantum output p x : 



W p x 



(1) 



Such a channel is known as a classical-quantum or cq channel 
for short. Two parameters that characterize the performance 
of this channel for transmitting classical data are the fidelity 



F (W) = Wy/Pny/PiWi an d the symmetric Holevo infor- 
mation I(W) = H((p Q + Pl )/2) - [H (p ) + H ( Pl )] /2, 
where H (a) = — Tr{ a log 2 er} is the von Neumann entropy. 
A channel is near perfect for transmitting classical data if 
I (W) ~ 1 and F (W) ~ 0, and it is nearly useless for this 
task if I (W) w and F (W) « 1 (Ref. [2| gives a precise 
relationship between these channel parameters). 

When encoding classical information for the channel W, 
we consider N = 2" uses of it, so that the resulting channel 
is of the form: x N = x\ ■ ■ ■ xjy — > p x N = p Xl ® • • • p XN , 
where x N is the classical input sequence and p x N is the output 
quantum state. Arikan's idea of channel combining and split- 
ting extends to this case by considering his encoding matrix 
Gn CD acting on an input sequence u N : u N —> p u n Gn , where 
u N Gm is the binary vector resulting from multiplication of the 
row vector u N by Arikan's encoding matrix Gn- We can then 

(i) 

define the "split" channels W N from the above "combined" 
channels as follows: 

W^:Ui^p^ B \ (2) 

i-1 
U l 

A ^Esm^v (4) 

The interpretation of this channel is that it is the one "seen" by 
the bit Ui if all of the previous bits u 1 ^ 1 are available and if we 
consider all the future bits uf +1 as randomized. This motivates 
the development of a quantum successive cancellation decoder 
J21 that attempts to distinguish Ui = from Uj = 1 by 
adaptively exploiting the results of previous measurements and 
quantum hypothesis tests for each bit decision. 

The above synthesized channels polarize, in the sense that 
some become near perfect for classical data transmission while 
the others become near useless. This result follows in the 
quantum case by a martingale argument similar to Arikan's: 
model the channel splitting and combining as a random 
birth process, and it is possible to prove that the parameters 
I(W$) and F(Wjy) converge almost surely to zero-one 
valued random variables in the limit of many recursions of 
the encoding. The following theorem characterizes the rate 
with which the channel polarization effect takes hold [8|, 0: 

Theorem 1: Given a binary input c q channel W and any 

(3 < 1/2, it holds that lim„^oo Pi i {\J F(W$) < 2- 2 "' 3 } = 
I (W), where n indicates the level of recursion for the 
encoding, W^} is a random variable characterizing the I th 
split channel, and F(W$) is the fidelity of that channel. 

Suppose now that we can determine which channels are 
the good ones and which ones are bad. In this case, we can 
construct a polar coding scheme by dividing the synthesized 
channels according to the following polar coding rule: 

Q N {W,P) e{ i£ [N] : Jf(W$) < 2-""}, (5) 

and Bjf(W,p) = [N] \ Gn(W,(3), so that G N (W,fi) is 
the set of "good" channels and Bn(W,/3) is the set of 



"bad" channels. The sender then transmits the information bits 
through the good channels and "frozen" bits through the bad 
ones. A helpful assumption for error analysis is that the frozen 
bits are chosen uniformly at random such that the sender 
and receiver both have access to these frozen bits. Ref. O 
provided an explicit construction of a quantum successive 
cancellation decoder that has an error probability equal to 
G (2-( 1 / 2 ) Ar ' 3 ) — let {aI^ 4 ^} denote the corresponding decod- 
ing positive operator-valued measure (POVM), with uj^ the 
information bits and u^c the frozen bits. 

II. Polar codes for quantum communication 

Wilde and Renes recently showed how to construct quantum 
polar codes for sending quantum data over an arbitrary qubit- 
input quantum channel [7| (before this, Renes et al. addressed 
the case of Pauli channels (6)). The polar coding scheme in 
Ref. [7] follows the general approach of Renes and Boileau 
(9), in which they showed how to build up quantum codes 
from classical codes for two different cq channels induced 
by the quantum channel of interest. We briefly overview the 
approach given in Ref. JTJ. 

Suppose that a quantum channel Af connects a sender to a 
receiver. Consider that the sender can induce a cq "amplitude- 
basis" channel Wa from M by inputting a computational basis 
state \z) depending on some classical bit z: 

W A :z^N{\z){z\). (6) 

Since the induced channel Wa is of the form in Q, it is 
possible to construct a polar code for it of rate / (WOi), 
consisting of the classical encoding Gm and a corresponding 
quantum successive cancellation decoder as described in the 
previous section. Now suppose that Alice shares a maximally 
entangled state |<£>) with Bob, where 

|$>=2- 1 /2 \*)\z). 

2£{0,1} 

Alice can modulate her share of this state with Z x depending 
on some bit x, so that the global state for Alice and Bob 
becomes 2^ 1 /' 2 (|0) |0) + (—l) x |1) |1)). If Alice sends her 
share of the state through the channel, this process induces 
a cq "phase-basis" channel Wp of the following form: 

Wp:x^{N®I){Z x \$>){>S>\Z x ). (7) 

It is also possible to construct a polar code for Wp of rate 
I (Wp), again consisting of the classical encoding Gn and a 
quantum successive cancellation decoder. 

The insight of Refs. H, Q is to build up a quantum polar 
code from the polar codes for the above two channels. The 
sender and receiver can compute offline which channel inputs 
will be good or bad for the amplitude and phase channels. 
The sender then inputs the following types of qubits into the 
different encoder inputs: 

1) Information qubits into the synthesized channels that are 
good in both amplitude and phase. 

2) Ancilla qubits in the state |0) into the synthesized 
channels that are bad in amplitude and good in phase. 



3) Ancilla qubits in the state |+) = (|0) + |1)) into 
the synthesized channels that are good in amplitude and 
bad in phase. 

4) Shares of ebits in the state |$) into the synthesized 
channels that are bad in both amplitude and phase. 

The intuition behind this approach is that the receiver should 
be able to recover quantum data if it is possible to recover 
classical data encoded into complementary variables. Thus, the 
sender puts information qubits into the synthesized channels 
which are good for both bases. The receiver will need help 
in the decoding process in the form of frozen ancilla qubits 
for the synthesized channels which are bad in one of the 
bases (this is the same as in Arikan's classical polar coding 
scheme fl], with the exception that the ancillas should be 
frozen in the basis that is bad). Finally, if both bases are 
bad, a shared ebit is in some sense frozen in both bases, 
because Bob can always predict the outcome of a Pauli X 
or Z measurement that Alice performs on her end if she tells 
him which measurement she performs. 

The net rate of this quantum polar coding scheme (rate 
of quantum communication minus the rate of entanglement 
consumption) is equal to the symmetric coherent information 
rate Q, @: 

I{W A ) + I(W P )-l = I{A)B) Mm 

= H ( B )aT(<S>) ~ H ( AB )n(1>) ■ 

The polar code operates as follows (see Ref. [7 1 for a detailed 
discussion of the operation). Alice sends qubits as described 
above into a coherent version of Arikan's encoder (this just 
consists of quantum CNOT gates). Bob then exploits a co- 
herent version of the quantum successive cancellation decoder 
for the amplitude cq channel in order to coherently decode the 
information qubits, placing them in a register C . Bob sends 
the coherently decoded information qubits through the inverse 
of the coherent polar encoder. This process induces the phase 
channel in (j7]) from the point of view of the phase basis of the 
qubits at the input of the encoder. Bob then acts with a co- 
herent version of the quantum successive cancellation decoder 
for the phase channels, coherently placing the outcomes in a 
register D . Finally, he performs a linear number of CNOT 
gates from the C N registers to the D N registers. The result is 
that he will have decoded the information qubits with a fidelity 
that is lower bounded by 1 — 2^/ea — 2-^/Fp, where €a and 
ep are the error probabilities of the respective polar codes for 
the cq amplitude and phase channels. 

The polar encoder is efficient because it is the same network 
of CNOT gates in Arikan's encoder (this part requires only 
O (N log N) operations). For general channels, it is still an 
open question to determine if the polar decoder has an efficient 
implementation. For Pauli channels or erasure channels, Renes 
et al. showed that the polar decoder has an efficient implemen- 
tation with O (NlogN) operations (6), essentially because the 
induced cq amplitude and phase channels are classical and they 
could thus exploit coherent versions of Arikan's efficient polar 
decoder in the quantum polar decoder. 



III. Polar codes for private communication 

An approach similar to the one above gives polar codes 
for sending classical data privately over a quantum wiretap 
channel Q. Again, this approach follows the general approach 
of Renes and Boileau [9| for transmitting private classical 
data, in which we build up such a protocol by considering two 
induced cq channels for complementary variables. A quantum 
wiretap channel j\f A ^ BE has an input system A for the sender, 
an output system B for the legitimate receiver, and an output 
system E for the wiretapper. The goal in this setting is for the 
sender to transmit classical data to the legitimate receiver in 
such a way that the wiretapper obtains a negligible amount of 
information about the input data. 

We now review the private polar coding protocol from 
Ref. 0. The first cq channel M.a that we consider has Alice 
prepare a quantum state p z depending upon a bit z and feed 
this state into the quantum wiretap channel Af: 

M A -z^M A ^ B { Pz ). (8) 

By purifying the input state p A to \ip z ) ASl and the quantum 
wiretap channel j\f A ^ BE to the isometric map jj A r > ' BES 2 
(where Si and S2 are "shield" systems not possessed by the 
wiretapper ifTOl . ifTTI '). we can see that the above channel 
arises by tracing over the wiretapper's system E and the shield 
systems S1S2: 

The channel in (jHJ is a cq channel, and as such, we can 
construct a polar code for it with rate I (Ma) along with 
an encoder and quantum successive cancellation decoder ||2) . 

The other cq channel that we consider is in some sense a 
virtual channel because we only make use of it in order to 
reason about the security of our private polar coding protocol 
(we do not actually exploit a decoder for it in the operation of 
our protocol). This channel is a phase channel with quantum 
side information. Suppose that Alice possesses an entangled 
state of the following form: 

\<p) CAS ^2-^Y,\*f\^) ASl - 

z 

Alice could then modulate the C system of the above state 
by applying the phase operator Z x to it, depending on some 
input bit x. Suppose then that she is able to transmit the A 
system through the channel M A ^ B to Bob and the C and Si 
systems through an identity channel. The resulting cq phase 
channel to Bob is as follows: 

M P : x -> M A ^ B ((Z X ) C \<p) (v\ CASl (Z x f). (9) 

Although we stated that this is a virtual channel, it is possible 
in principle to construct a polar code for it with rate I(A4p) 
along with a quantum successive cancellation decoder Q. 

The usefulness of the phase channels for privacy may not 
be readily apparent, but an uncertainty relation from Ref. |9| 
clarifies the link. Indeed, consider the following channel to the 
wiretapper: 

M E -z^N A ^ E ( Pz ). (10) 



The uncertainty relation (Lemma 2 of Ref. [9]) for our case 
is as follows: 

I(M p ) + I(Me) = 1. (11) 

The implication of this uncertainty relation is that if the phase 
channel Mp to Bob is nearly perfect, then the amplitude 
channel Me is nearly useless to the wiretapper and vice versa. 

We can then exploit the above uncertainty relation to con- 
struct a polar code for private communication. The structure 
is similar to that in the previous section, though the resources 
used are in some sense classical versions of the quantum 
resources. The sender and receiver compute offline which 
channel inputs will be good or bad for the amplitude or phase 
channels in ([8]) and (|9j, respectively. The sender then inputs 
the following types of qubits into the different encoder inputs: 

1) Information bits to be kept private into the synthesized 
channels that are good in both amplitude and phase. 

2) Ancilla bits initialized to into the synthesized channels 
that are bad in amplitude and good in phase. 

3) Ancilla bits that are randomized (0 or 1 with probability 
1/2) into the synthesized channels that are good in 
amplitude and bad in phase. 

4) Shares of secret key bits (0 or 1 with probability 1/2 
and known to both Alice and Bob) into the synthesized 
channels that are bad in both amplitude and phase. 

The security of this approach follows from the uncertainty 
relation in (111. If Alice places the information bits into 



and outputs this state with probability 1 — e or an orthogonal 
erasure flag |e) with the complementary probability: 



the synthesized channels that are good in amplitude and 
phase, then Bob will be able to recover them reliably since 
the channels are good in amplitude and Eve learns only 
a negligible amount of information about them due to the 



uncertainty relation in (111 and the fact that the channels are 
good in phase for Bob. The sender places ancilla bits set to 
in the channels that are bad in amplitude and good in phase in 
order to help Bob in decoding the information bits. She places 
randomized bits into the channels that are good in amplitude 
and bad in phase in order randomize Eve's knowledge of them. 
Finally, she places secret key bits into the channels that are 
bad for both variables in order to help Bob decode and to 
randomize Eve's knowledge of the information bits. 

The net rate of this private polar coding scheme (rate of 
private classical communication minus the rate of secret key 
consumption) is equal to the symmetric private information 
rate 0, 0: 

/ (M A ) + I (M P ) -1 = 1 (M A ) - I (M E ) ■ 

We can estimate the reliability and security of this scheme by 
considering the performance of the constituent cq polar codes. 

IV. Rigorous Bounds on Performance for the 
Quantum Erasure Channel 

We now arrive at our first result, where we provide rigorous 
bounds on the performance of quantum polar codes on the 
quantum erasure channel. Our development here is both simi- 
lar to and extends that in Section V-D of Arikan [ 1 1. Recall that 
a quantum erasure channel takes a qubit in the state p as input 



p^(l-e)p + e\e) (e| 



(12) 



First, consider that in any given "run" of a quantum polar 
code of the form reviewed in Section [IIJ the block error 
probability P e is as follows: 

P e = Pr {a en } + Pr {0 OT } - Pr {a en , <f> m } , (13) 

where a elT and </> elT are the events that an error occurs during 
the decoding of the induced amplitude and phase channel, 
respectively. (The above formula ensures that we count a 
quantum block error when either an amplitude or phase error 
occurs, but that we do not overcount when both errors occur.) 
By exploiting the law of total probability (that Pr{</> elT } = 
Pr {a elT , 0err} + Pr {<3ern 0err}> where a^- is the event that 
there is no amplitude decoding error), we have that 



P e = Pr{cw} + Pr{ot 

= Pr Kn-} + Pr {</>e, 
= Pr{a en -}+Pr{0 e , 

= Pr Kn-} + Pr {</>« 



tea} 

a^} Pr {a^-} 

a^} [1 - Pr{a elT }] 

a^} - Pr {0 err | a^} Pr {a en } . 



When considering performance for the quantum erasure 
channel, there is a symmetry that is helpful in simplifying the 
above expression. Consider that both the induced amplitude 



and phase channels for the quantum erasure channel in (12i 
are classical erasure channels with erasure probability e. Let 
Wjy be the i lh synthesized channel from the classical erasure 
channel, r\ a threshold parameter between and 1, and A ij]) 
the set of information bits chosen by Arikan's polar coding 
rule: A(r\) = {i : Z(W$) < t]}, with Z the Bhattacharya 
parameter. Recall that the order of these indices i is preserved 
for the induced amplitude channels but reversed for the in- 
duced phase channels. If we assume that r\ is the same for both 
the induced amplitude and phase channels, then it follows that 

Pr{a err } = Pr{0 elT | a^-} . 

This assumption is reasonable because both of the induced 
channels are erasure channels with the same erasure probabil- 
ity. By exploiting the above symmetry, we have that 

P e = Pr{a err }(2-Pr{a e rr}). 

Thus, in order to bound the block error probability P e 
from both above and below, we only need to bound Pr {a elT } 
from above and below. In order to do so, we can exploit the 
following upper bound from Section V-B of Arikan £[]: 



Pr{aerr} < U fr) 



E 

iGA(r]) 



z(w£>). 



Then by using Arikan's recursive equations for the erasure 
channel in (38) of Ref. JT], we can easily compute this upper 
bound on Pr{a err }. 

For a rigorous lower bound on Pr{a err }, we appeal to 
arguments similar to those in Section V of Ref. |2j, since 



one can always embed classical systems into quantum systems 
(also, recall that Arikan only provided a heuristic lower bound 
in Section V-D of Ref. [T)). We can write the error probability 
Pr{a err } as follows: 



-y 



2 N 



1 -Tr 



V), 



■ n 



(i),u\ 



• TP 3 

n (l),«l 



Pu N 



n 



• n 



(i),u\ 



n 



f -\t\ N—l 
{N),u 1 u N 



(14) 



where p u N is the encoded state and 11^, 



tB'' 



\ B " n 

are projectors corresponding to quantum hypothesis tests to 
decode each bit m. In the classical case, p u N is just a 
distribution given by the encoding and channel and EK-v Ui , 

H?„„ are indicator functions corresponding to 

{N),u{ un 

likelihood ratio tests. Consider the following operator inequal- 
ities which hold for commuting operators P x and P 2 such that 

0<Pi,P 2 <I: 

I-P 1 P 2 P 1 > J-Pi, 
/ - P1P2P1 > I - P 2 . 
We can exploit these recursively in order to obtain the follow- 



ing lower bound on the error term in ( [14) , since the projectors 

t b"' 



ttB" 

n (i),«i' ' 
channel: 



■ , nf Ar , all commute for a classical erasure 



H 2 i-i XI 2 Tr 1 ( 7 ^Mr 1 ".) X 2 N 

f/ r/'-iRN\ rri-lpfl 



= Pen- (Ui) 



where in the last line we have exploited the states and notation 
defined in Section V of Ref. Q. By recognizing that the 
quantity above is equivalent to the error in discriminating the 
states p,i Q and p. I 1 in a quantum hypothesis test and 
recalling the relationship between hypothesis testing error, the 
trace distance, and the fidelity (Bhattacharya parameter), we 
have the following lower bound from Ref. |12|: 

p e rA^)>\(i-\/i-z(wtfy 

Thus, since this bound holds for any index i, we obtain the 
following lower bound on Pr {a err } for the quantum erasure 
channel: 

Pr {a elr } > L(rj) = max - [ \ - 



1 - z{w$y 



i€A(v) 2 

This bound is also easy to compute by exploiting Arikan's 
recursive equations for the erasure channel in (38) of Ref. (T). 

The above development then leads us to the following 
theorem: 

Theorem 2: The block error probability P e for a quantum 
polar code of the form reviewed in Section [TT] when used on 
a quantum erasure channel is bounded as follows: 

L( V ) (2 - U ( V )) < P e < U ( V ) (2-1,(7?)). 
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Fig. 1 . The figure depicts the upper (blue solid curves) and lower (red dashed 
curves) bounds from Theorem [2] for a quantum erasure channel with erasure 
probability e = 0.15. These bounds determine the trade-off between block 
error probability and rate for quantum polar codes designed for a quantum 
erasure channel. The channel's quantum capacity is equal to 1 — 2e = 0.7 
qubits per channel use. Next to each curve, we also show the block length 
for a quantum polar code to which the corresponding bound applies. 



Figure [T] depicts the bounds given by Theorem [2] What 
we observe is that a large block length « 2 20 is required in 
order for a quantum polar code to operate in a regime near the 
quantum capacity. Also, the margins between the upper and 
lower bounds from Theorem|2]are rather large. Thus, it is clear 
that numerical simulations would be helpful in determining 
more accurate estimates of the performance of quantum polar 
codes on the quantum erasure channel. 

V. Numerical simulation results 

We now describe the results of our numerical simulations 
of quantum polar codes on quantum erasure, depolarizing, and 
BB84 channels. Our simulations rely on two observations from 
Ref. (6] (and noted in the previous section): 

• The cq channels in |6} and |7]i corresponding to the 
respective induced phase and amplitude channels result 
in distinguishable states at the output (for erasure, depo- 
larizing, and BB84 channels). 

• The performance of the resulting quantum polar code is 
directly related to the performance of polar codes for the 
constituent amplitude and phase channels. 

The first observation implies that the induced channels 
are effectively classical. Specifically, for the quantum erasure 
channel with erasure probability e, both the induced amplitude 
and phase channel correspond to a classical erasure channel 
with erasure probability e. The second observation implies that 
we can simulate the performance of a quantum polar code for 
Pauli or erasure channels by simulating the performance of 
Arikan's successive cancellation decoder for the constituent 
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Fig. 2. Diagram of the information bits for the induced amplitude and phase 
channels for code length N = 8 and a classical rate Rq = 0.75. The numbers 
represent the symmetric capacities for each virtual channel when e = 0.15. 
The blue boxes represent the 'good' channels, which for our example are 
the virtual channels with symmetric capacities above the lowest quartile. The 
quantum code then chooses the four bits which are good for both amplitude 
and phase, resulting in a rate Rq = 0.5 quantum polar code. 



channels. Importantly, this simulation is efficient because 
Arikan's successive cancellation decoder is efficient, requiring 
only 0(N log N) operations. 

Our Monte Carlo simulation proceeds by randomly gener- 
ating information bits, encoding them with the polar encoder, 
transmitting them over a given channel, and then decoding 
with the successive cancellation decoder. In our simulations, 
we vary the choice of the classical rates and R^ mp ^ 

for the respective phase and amplitude channels. For the case 
of the quantum erasure channel, it is optimal to choose these 
rates to be equal because the two induced amplitude and phase 
channels are classical erasure channels with the same erasure 
probability and hence have the same capacity. 

Figure [2] depicts the process by which we choose an 
example quantum polar code with N — 8 and classical rates 



R 



(amp) 



= R 



(Ph) 



= 0.75. The virtual channels' symmetric 



capacities in the erasure case can be calculated explicitly 
and efficiently via the recursive formulas in (38) of Ref. (T). 
The reason that the symmetric capacities are reversed for the 
induced phase channel is that the action of the quantum CNOT 
gates in the encoding are reversed when acting on the phase 
basis (as first observed in O). We end up with four channels 
"good" for both amplitude and phase and Rq = 4/8 = 0.5. 
In the case where the bad virtual channels are completely 
clustered at the ends, it holds that Rq = R { g b) + R^ mp) - 1. 
However, we note that the symmetric capacities of the virtual 
channels are generally not ordered sequentially {e.g., u$ has a 
larger symmetric capacity than 114) and in larger block lengths 
N, the good virtual channels become distributed more non- 
contiguously, resulting in slightly higher quantum rates. 

After randomly choosing the information bits u^, the 
encoding is performed by multiplying by the encoding matrix 
Gn DD> and the codeword is transmitted through the channel 
according to its probability transition matrix. For the amplitude 
channels, we then sequentially decode u\, ita, • • • ■ We always 
decode any frozen bits correctly, and we decode the informa- 
tion bits to be either '0' or '1' according to the likelihood 
ratio, which for each bit is calculated recursively in 0(log N) 
steps by incorporating the values of previously decoded bits 
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Fig. 3. Block error rate P e versus quantum data rate for a quantum erasure 
channel with e = 0.15 and different block lengths. Yellow, red and blue 
points show, respectively, N = 64, 256, and 1024. Points show simulated 
results and lines are linear interpolations between points. The error bars here 
and in subsequent figures reflect the 90% confidence interval for the error 
rates given the number of simulations performed and errors observed. The 
quantum capacity of this channel is equal to 0.7 qubits per channel use. 



HI. For the induced phase channels, we decode in the reverse 
order, due to the reversal mentioned above. 

We calculated the block error rates Pr {a err } and 
Pr {<f>err | ct en } for the respective induced amplitude and phase 
channels separately by performing M such simulations of 
each, in each case comparing with the actual encoded word. 
The block error rate for quantum data transmission is equal 



to ( 1 3 1 by counting an error whenever either the amplitude or 
phase decoding fails (but not overcounting when both fail). 
For most cases, we ran M — 50, 000 simulations in order to 
get statistically significant results down to P e < 10~ 4 (though 
for some larger P e cases we ran M = 5, 000, as this was more 
than sufficient to obtain small statisical error bars). 

A. Results for Quantum Erasure Channels 

Figure [3] depicts the block error rate versus quantum data 
rate for three different code lengths N — 64, 256, and 1024. 
All of the block error rates start to bend down slightly below 
the quantum capacity for the channel (which in this case is 
equal to 0.7 qubits per channel use), and the longer block 
lengths improve more rapidly at lower rates. Figure |4ja) then 
plots the block error rate versus erasure probability e for one 
particular quantum data rate Rq w 0.4. Note that the exact 
rates for different block lengths N are slightly different due to 
a slightly increasing fraction of overlapping phase-good and 
amplitude-good channels at larger block lengths. 

Figure |4|b) summarizes performance for a variety of erasure 
probabilities e and quantum data rates Rq by plotting the 
threshold rate at which P e < 10~ 4 . It also plots the quantum 
capacity 1 — 2e of the quantum erasure channel [13|. We 
observe that for the selected block lengths, quantum polar 
codes can perform increasingly near to capacity as the erasure 
probability e decreases. 
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Fig. 4. Quantum erasure channel simulation results, (a) Block error rate 
P e versus erasure probability e for N = 64 (yellow), N = 256 (red), and 
N = 1024 (blue) and quantum rate Rq = 0.398 (the N = 64 points were 
for Rq = 0.375). (b) Threshold quantum data rate Rq for block error rate 
P e < 10 — 4 versus erasure probability t. The black curve shows the quantum 
capacity 1 — 2e. 



B. Results for Quantum Depolarizing Channels 

We also simulated the performance of quantum polar codes 
for the quantum depolarizing channel with depolarizing prob- 
ability p, modeled as follows: 

p^ (l-p)p + p/3 (XpX + YpY + ZpZ). 

In this case, the induced amplitude channel has orthogonal 
outputs, with errors occuring only in the case of an X or 
Y flip. Thus, it is equivalent to a classical binary symmetric 
channel (BSC) with flip probability 2p/3. The induced phase 
channel has four possible outputs (corresponding to the four 
orthogonal Bell states), and it is thus equivalent to a classical 
channel with the following probability transition matrix: 

1-p p/3 p/3 p/3 
p/3 p/3 p/3 1-p 

One important subtlety for both the induced amplitude and 
phase channel for the quantum depolarizing channel is that 
there is no explicit method for calculating the symmetric 
capacities of each virtual channel (as in the case of an erasure 
channel). From Ref. (TJ, we know that there is a choice of 
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Fig. 5. Quantum depolarizing channel simulation results, (a) Block error rate 
P e versus physical error probability p for N = 64 (yellow), N = 256 (red), 
and N = 1024 (blue) and quantum data rate Rq = 0.30. (b) Threshold 
Rq for P e < 10 -4 versus p. The black curve shows the hashing limit 
l-tf 2 (p)-plog 2 (3). 

good channels with exponentially decreasing error probability 
for any rate below capacity, but Ref. Qj does not give an 
explicit prescription for choosing them. One strategy from 
Ref. [ 14], which we employ here, is to calculate the symmetric 
capacity S of a channel and use an effective erasure probability 
e = 1 — S to choose the good channels according the 
prescription for erasure channels. It is not fully known how 
close to optimal this heuristic is. 

Another important subtlety is that the induced amplitude 
and phase channels have different symmetric capacities (with 
the phase channel's being slightly larger). For this reason, it 
is better to choose the amplitude classical rate slightly lower 

To account for this, we 



than the phase rate: R ( * mp > < R^ h) . 
varied the rates R^, and R^, ' independently and chose 
the combination of these two rates which optimized P e for 
a given quantum data rate Rq. We found empirically that 



choosing R^ m? ^ 



gave near optimal results across 
the channel parameters and rates we considered but varying 
around this choice gives slight improvements. 

Figure |5ja) plots the block error rate P e for polar codes 
with quantum data rate Rq « 0.30, and Figure |3Jb) plots 
the threshold for which P e < 10~ 4 along with the hashing 
bound 1 — H2 (p) — p log 2 (3) (this is the rate achievable by a 
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Fig. 6. BB84 channel simulation results. Threshold Rq for P e < 10 4 
versus flip probability /. The black curve is the hashing limit 1 — 2H2(f). 



random stabilizer code on the depolarizing channel). We see 
the performance relative to the bounds is not as strong as for 
the quantum erasure channel. Nevertheless, we obtain fairly 
large quantum data rates with low block error rate P e and 
moderate block lengths TV = 1024. 

In Ref. 11151 . Kasai et al. presented results for quantum codes 
with code lengths in the range N — 4000-35, 406 and achieved 
quantum rates higher than ours at the P e < 10~ 4 threshold. We 
should perform simulations of higher N to compare our coding 
scheme with theirs. However, it appears from our evidence 
so far that at these code lengths, achievable quantum polar 
code rates will be lower than theirs, unless we exploit some 
improved method such as that in Ref. [16] for computing the 
good and bad virtual channels. 

C. Results for the BB84 Channel 

Finally, we consider performance on the "BB84 channel" 
[17 1, which is a concatenation of a bit-flip and phase-flip 
channel, each having some equal flip probability / (see 
Section 2.1.1 of Ref. ifTTIl for the relevance of this channel in 
quantum key distribution). For the BB84 channel, the induced 
amplitude and phase channels are simply binary symmetric 
channels, each with error probability /. Because the two 
channels are equivalent, as in the erasure channel case, it is 



optimal to choose equal underlying rates: R 



(amp) 
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(phase) 



Figure [6] plots our results for quantum data rates achieving the 
P e < 10~ 4 threshold. We see quantum polar code performance 
is very similar to that for the depolarizing channel. 

VI. Application to Private Communication 

Our simulations above serve equally well as simulations of 
polar codes for private classical data transmission over the 
channels we considered (where here the eavesdropper gets 
access to everything that the receiver does not obtain). This 
follows essentially from Theorem 1 of Renes fl8l . in which he 
proves that a lower bound on Bob's guessing probability for 
the phase variable serves as an upper bound for security against 
Eve with respect to the amplitude variable. Thus, the parameter 



the protection of classical data sent through the channels we 
have considered. 

VII. Discussion 

In summary, we have bounded the performance and con- 
ducted simulations of quantum communication over the quan- 
tum erasure channel, by utilizing the recently proposed quan- 
tum polar codes in Refs. 0, Q. We found high quantum 
data rates for moderate block lengths N = 1024. We also 
performed simulations of these codes over the depolarizing 
and BB84 channels. In these cases, we found that quantum 
polar codes for these block lengths still performed ably but 
gave performance somewhat further from the hashing limits. 

Going forward from here, it is important to explore the 
performance of larger block lengths in order to compare 
the performance of quantum polar codes with other error 
correction schemes. In the depolarizing and BB84 cases, it 
would be worthwhile to carry out an analysis of the optimality 
of the erasure-matched code choice. Also, one could exploit 
the techniques from Ref. Ifl6l in order to choose which virtual 
channels to send the information bits through. 
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